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ABSTRACT 

We report on our serendipitous pre-discovery detection and follow-up observations of the broad-lined Type Ic 
supernova (SN Ic) 2010ay at z — 0.067 imaged by the Pan-STARRSl 37r survey just ~4 days after explosion. The 
supernova (SN) had a peak luminosity. Mg ~ —20.2 mag, significantly more luminous than known GRB-SNe and 
one of the most luminous SNe Ib/c ever discovered. The absorption velocity of SN 2010ay is usi ~ 1 9 x 10 3 km s _1 
at ~40 days after explosion, 2-5 times higher than other broad-lined SNe and similar to the GRB-SN 2010bh 
at comparable epochs. Moreover, the velocity declines ~2 times slower than other SNe Ic-BL and GRB-SNe. 
Assuming that the optical emission is powered by radioactive decay, the peak magnitude implies the synthesis 
of an unusually large mass of 56 Ni, My, = 0.9 M 0 . Applying scaling relations to the light curve, we estimate a 
total ejecta mass, M e j 4.7 M 0 , and total kinetic energy, E K ~ 1 I x 10 51 erg. The ratio of M N ; to M c] is ~2 
times as large for SN 2010ay as typical GRB-SNe and may suggest an additional energy reservoir. The metallicity 
(Iog(0/H)ppo4 + 12 = 8.19) of the explosion site within the host galaxy places SN 2010ay in the low-metallicity 
regime populated by GRB-SNe, and ~0.5(0.2) dex lower than that typically measured for the host environments 
of normal (broad-lined) SNe Ic. We constrain any gamma-ray emission with Ey< 6 X 10 48 erg (25-150 keV), 
and our deep radio follow-up observations with the Expanded Very Large Array rule out relativistic ejecta with 
energy E > 10 48 erg. We therefore rule out the association of a relativistic outflow like those that accompanied 
SN 1998bw and traditional long-duration gamma-ray bursts (GRBs), but we place less-stringent constraints on a 
weak afterglow like that seen from XRF 060218. If this SN did not harbor a GRB, these observations challenge the 
importance of progenitor metallicity for the production of relativistic ejecta and suggest that other parameters also 
play a key role. 
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1. INTRODUCTION 

Recent observations have shown that long-duration gamma- 
ray bursts (GRBs) are accompanied by Type Ic supernovae 
(SNe Ic) with broad absorption features (hereafter “broad- 
lined,” BL) indicative of high photospheric expansion velocities 
(see Woosley & Bloom 2006, for a review). This GRB-SN 


18 L.C. is a Jansky Fellow of the National Radio Astronomy Observatory. 

19 Emeritus 


connection is popularly explained by the favored “collapsar 
model” (MacFadyen et al. 2001) in which the gravitational 
collapse of a massive (M > 20 M Q ) progenitor star gives birth 
to a central engine — a rapidly rotating and accreting compact 
object — that powers a relativistic outflow. At the same time, 
not all SNe Ic-BL show evidence for a central engine. Radio 
observations constrain the fraction of SNe Ic-BL harboring 
relativistic outflows to be less than a third (Soderberg et al. 
2006a, 2010). 
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The physical parameter(s) that distinguishes the progenitors 
of GRB-associated SNe from other SNe Ic-BL remains debated, 
while theoretical considerations indicate that progenitor metal- 
licity may play a key role (Woosley & Heger 2006). In the col- 
lapsar scenario, massive progenitor stars with metallicity above 
a threshold Z > 0.3 Z Q lose angular momentum to metal-line- 
driven winds, preventing the formation of a rapidly rotating 
compact remnant and, in turn, a relativistic outflow. At the same 
time, the hydrogen-free spectra of SNe Ic-BL indicate that their 
stellar envelopes have been stripped prior to explosion, requir- 
ing higher metallicities (e.g., Z « Z 0 ) if due to radiation- 
driven winds (Woosley et al. 1995). Alternatively, short-period 
(~0.1 days) binary interaction may be invoked to spin up stars 
via tidal forces, as well as cause mass loss via Roche lobe over- 
flow (Podsiadlowski et al. 2004; Fryer & Heger 2005). However, 
even in the binary scenario, GRB formation is predicted to occur 
at higher rates in lower metallicity environments, where the ra- 
dius and mass-loss rates of stars should be smaller (Izzard et al. 
2004). 

Observationally, most GRB-SNe are discovered within dwarf 
star-forming galaxies (Fruchter et al. 2006) characterized by 
sub-solar metallicities, Z < 0.5 Z G (Levesque et al. 2010a). 
This has been interpreted as observational support for the 
metallicity-dependent collapsar model (Stanek et al. 2006). 
Meanwhile, it was noted early on that SNe Ic-BL without asso- 
ciated GRBs have been found in more enriched environments 
(e.g., Modjaz et al. 2008; Prieto et al. 2008), motivating the 
suggestion of an observationally determined cutoff metallicity 
above which GRB-SNe do not form (Kocevski et al. 2009). 
This metallicity difference may be partly attributable to the 
different survey strategies: SNe have been primarily discov- 
ered by galaxy-targeted surveys biased toward more luminous 
(and therefore higher metallicity) environments, while GRB 
host galaxies are found in an untargeted manner through their 
gamma-ray emission (Stanek et al. 2006). 

Against this backdrop of progress, recent observations have 
begun to call into question some aspects of this scenario. 
First, there are two reported long-duration GRBs in solar or 
super-solar metallicity environments (GRB 020819, Levesque 
et al. 2010c; GRB 051022, Graham et al. 2009). Similarly, the 
luminous radio emission seen from SN Ic-BL 2009bb pointed 
unequivocally to the production of copious relativistic ejecta 
resembling a GRB afterglow (Soderberg et al. 2010), while 
the explosion environment was characterized by a super-solar 
metallicity, Z ~ 1-2 Z Q (Levesque et al. 2010b). Together 
with the growing lack of evidence for massive progenitor 
stars for SNe Ic in pre-explosion Hubble Space Telescope 
images (Smartt 2009), a lower mass (M ~ 10-20 M Q ) binary 
progenitor system model (with a gentler metallicity dependence) 
is gaining increasing popularity (Yoon et al. 2010). Multi- 
wavelength studies of SNe Ic-BL in metal-poor environments 
may shed further light on the role of metallicity in the nature 
of the progenitor and the explosion properties, including the 
production of relativistic ejecta. 

Fortunately, with the recent advent of wide-held optical 
transient surveys (e.g., Catalina Real-time Transient Survey 
(CRTS), Drake et al. 2009; Panchromatic Survey Telescope 
and Rapid Response System (Pan-STARRS 1 , abbreviated PS 1 ), 
Kaiser et al. 2002; Palomar Transient Factory, Law et al. 
2009) SNe Ic-BL are now being discovered in metal-poor 
environments, Z ~ 0.5 Z Q , thanks to an unbiased search 
technique (see, e.g., Stanek et al. 2006; Modjaz et al. 2008, 2011; 
Young et al. 2008; Arcavi et al. 2010; Leloudas et al. 201 1; Kelly 


Table 1 

SN 2010ay Host Galaxy SDSS J123527.19+270402.7 


Parameter 


Value 

R.A. 


12 h 35 m 27 s 19 02000) 

Deck 


+27°04'02'.'7 (J2000) 

Redshift (z) 


0.0671 

Petrosian radius 


17355 

Photometry 2 

u! 


19.56 ± 0.03 mag 

&' 


19.02 ± 0.01 mag 

r' 


19.02 ± 0.01 mag 

i' 


18.69 + 0.01 mag 

z' 


18.87 ±0.04 mag 

U 


19.50 ± 0.06 mag 

B 


19.02 ± 0.05 mag 

V 


19.13 ±0.05 mag 

R 


18.94 + 0.05 mag 

I 


18.90 ± 0.06 mag 

Extinction 

E(B - V)MW b 


0.017 mag 

E(B - R)h ost c 


0.2 mag 


Notes. SDSS host galaxy properties and ugriz photometry. 
a Model magnitudes from SDSS DR8 (Aihara et al. 2011). 

Host galaxy photometry has not been corrected for extinction. 

UBVRI photometry has been converted from the SDSS ugriz 
measurements using the transformation of Blanton & Roweis 
(2007). 

b The Milky Way extinction as determined by Schlegel et al. 

(1998), assuming Ry =3.1. 

c The host galaxy extinction determined from the SDSS 
spectrum centered on the galaxy nucleus, via the Balmer 
decrement as described in Section 2.3. 1 . 

& Kirshner 2011; Sanders et al. 2012; Stoll et al. 2012). In this 
paper, we present pre-discovery Pan-STARRS 1 imaging and 
multi-wavelength follow-up observations for the broad-lined 
Type Ic SN 2010ay discovered by CRTS (Drake et al. 2010). In 
Section 2, we report our optical (Pan-STARRS 1 , Gemini, and 
William Herschel Telescope (WHT)) and radio (Expanded Very 
Large Array (EVLA)) observations. In Section 3, we model the 
optical light curve and analyze the spectra to derive the explosion 
properties of SN 20 1 Oay. In Section 4, we use our observations of 
SN 201 Oay with the EVLA to place strict limits on the presence 
of relativistic outflow. In Section 5, we draw from gamma-ray 
satellite coverage to rule out a detected GRB in association 
with SN 201 Oay. In Section 6, we derive the explosion site 
metallicity and find it to be significantly sub-solar and typical 
of most GRB-SN host environments. In Section 7, we discuss 
the implications of these findings in the context of the explosion 
and progenitor properties and conclude in Section 8. 

2. OBSERVATIONS 
2.1. Discovery by CRTS 

SN 201 Oay was discovered by the CRTS (Drake et al. 2009) 
on 2010 March 17.38 UT (Drake et al. 2010) and designated 
CSS100317:123527+270403, with an unfiltered magnitude of 
m ~ 17.5 mag and located <1" of the center of a compact 
galaxy, SDSS J123527. 19+270402.7, at z = 0.067 (Table 1). 
We adopt a distance Di = 297.9 Mpc to the host galaxy 20 


20 We assume Ho = 71 km s 1 Mpc 1 , Q\ = 0.73, and Q.m = 0.27 . 
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Figure 1. Images illustrating the Pan-STARRSl pre-discovery detection of SN 2010ay and the surrounding field, (a) Pre-explosion gr /-composite image from the 
SDSS DR7 (Abazajian et al. 2009), observed 2004 December 21, (b) SDSS /-band image geometrically registered to the PS1 image frames and zoomed in to the host 
galaxy of SN 2010ay, (c) /pi-band image from the 3 tt survey of PS1, observed 2010 February 25, (d) the difference of the SDSS /' and PS1 /pi images. The transient 
emission can be seen in the NE corner of the last panel. Nearby stars are included in the zoomed-in frames to illustrate the efficiency of subtraction. 

(A color version of this figure is available in the online journal.) 


and note that the Galactic extinction toward this galaxy is 
E(B — V) = 0.017 (Schlegel et al. 1998). Pre-discovery 
unfiltered images from CRTS revealed an earlier detection of 
the SN on Mar 5.45 UT at m ~ 18.2 mag and a non-detection 
from February 17.45 UT at m > 18.3 mag (Drake et al. 2010). 

A spectrum obtained on March 22 UT revealed the SN 
to be of Type Ic with broad features, similar to the GRB- 
associated SN 1998bw spectrum obtained near maximum light 
(Filippenko et al. 2010). This classification was confirmed 
by Prieto (2010), who additionally reported photometry for 
the SN (see Table 2). After numerically subtracting the host 
galaxy emission, they estimate an unusually luminous absolute 
magnitude of V & —19.4 mag. 

2.2. Pre-discovery Detection with Pan-STARRSl 3tt 

The field of SN 2010ay was serendipitously observed with 
the PS1 3ic survey in the weeks preceding its discovery. Pan- 
STARRS lisa wide-held imaging system at Haleakala, Hawaii, 
dedicated to survey observations (Kaiser et al. 2002). The PS1 
optical design (Hodapp et al. 2004) uses a 1.8 m diameter //4.4 
primary mirror and a 0.9 m secondary. The telescope illuminates 
a diameter of 3? 3. The Pan-STARRSl imager (Tonry & Onaka 
2009) comprises a total of 60 4800 x 4800 pixel detectors, with 
10 jum pixels that subtend 0.258 arcsec. The PS1 observations 
are obtained through a set of five broadband filters designated as 
gPi, rpi, ip | , zpi, and ypi. These filters are similar to those used 
in previous surveys, such as Sloan Digital Sky Survey (SDSS; 
Fukugita et al. 1996). However, the g pi filter extends 20 nm 
redward of g', the z pi filter is cutoff at 930 nm, and SDSS has 
no corresponding vpi filter (Stubbs et al. 2010). 

The held of SN 2010ay was observed on 2010 February 21 
(rpi band) and February 25 (ipi band; Figure 1). On each night, 
four exposures were collected following the strategy of the PS 1- 
3 jt survey (K. C. Chambers et al. 2012, in preparation). Follow- 
ing the CRTS discovery and announcement of SN 2010ay, we 
geometrically registered SDSS pre-explosion images to the PS1 


Table 2 

SN 2010ay Light-curve Photometry 


UT Date 

MJD 

/peak 3 

Filter 

m b 

m r c 

Source 

733820.4 

55244.4 

-29 


< 18.3 

<- 19.2 ± 0.2 

d 

733824.6 

55248.6 

-25 

rp\ 

< 22.0 ± 0.1 

<- 16.0 ± 0.1 

e 

733828.2 

55252.2 

-21 

*Pl 

21.1 ± 0.3 

- 16.8 ± 0.3 

e 

733836.4 

55260.4 

-13 


18.2 

- 19.3 ± 0.2 

d 

733848.4 

55272.4 

-1 


17.5 

- 20.4 ± 0.2 

d 

733853.2 

55277.2 

4 

B 

18.39 ± 0.05 


f 

733853.2 

55277.2 

4 

V 

17.61 ± 0.05 


f 

733853.2 

55277.2 

4 

R 

17.44 ± 0.05 

- 20.22 ± 0.07 

f 

733863.0 

55287.0 

14 

R 

18.2 ± 0.2 

- 19.2 ± 0.2 

g 

733873.0 

55297.0 

24 

g 

18.9 ± 0.1 


h 

733873.0 

55297.0 

24 

r 

18.3 ± 0.1 

- 19.0 ± 0.1 

h 

733873.0 

55297.0 

24 

i 

18.0 ± 0.1 


h 

734221.0 

55645.0 

372 

*Pl 

< 22.2 ± 0.2 

<- 16.22 ± 0.10 

e 

734225.0 

55649.0 

376 

rpi 

< 21.9 ± 0.1 

<- 15.7 ± 0.1 

e 


Notes. 

a Time since peak in days, relative to the fitted value: 2010 March 18.5 ± 0.2. 
b The measured apparent magnitude of the source, in the filter noted and without 
extinction correction. For the Pan-STARRS 1 photometry, a template image was 
subtracted; for the other points, the host galaxy flux has not been subtracted. 
c The absolute R magnitude of the SN. Filter conversion, host flux subtraction, 
and extinction correction have been performed as described in Section 3.1. 
d Drake et al. (2010), unfiltered (synthetic V band). 
e Photometry from the Pan-STARRS 1 3ic survey. 
f Prieto (2010). 

g Synthetic photometry obtained from our WHT spectrum as described in 
Section 2.3. 

h Photometry from our Gemini/GMOS observations described in Section 2.2. 


images and performed digital image subtraction using the ISIS 
package (Alard 2000). Photometric calibration was performed 
using many field stars in the SDSS catalog with magnitudes 
r « 13-19 mag. No residual flux was found in the difference 
r P1 -band image from 2010 February 21 with an upper limit of 
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Figure 2. Optical /?-band light curve of SN 2010ay, as compiled in Table 2 from 
CBET 2224 (red squares), the PS1 survey (blue), our Gemini images (black), 

and synthetic photometry (Section 2.3) based on our WHT spectrum (gold). 
Triangles denote upper limits. The thick dashed line represents the luminosity 
of an expanding fireball fit to our early-time photometry (Section 3.1). The thin 
dashed line is the SN Ib/c light-curve template of Drout et al. (201 1), and the 
gray held represents the standard deviation among its constituent photometry. 
The template is stretched by (1 + z) = 1.067 with the best-fit parameters 
hipeak = March 18 ± 2 (MJD 55273 d= 2) and MR pea k = —20.2 d= 0.2 mag. 
The solid line is the light curve of the GRB-SN 1998bw (Galama et al. 1998), 
redshifted to match SN 2010ay. The vertical dotted lines mark the epochs of our 
Gemini and WHT spectroscopy. 

(A color version of this figure is available in the online journal.) 

r > 22.0 mag. However, we detect residual flux at the posi- 
tion of SN 2010ay in the 4. | -band residual image from 2010 
February 25 with a magnitude of i' = 21.1 ± 0.3 mag. 

The field was again observed in the z P | band on 2011 March 
25 and the r Pi band on 2011 March 29, but no residual flux 
was detected in the subtractions at the SN position to limits of 
i' > 22.2 and r' > 21.9. 

In Table 2 and Figure 2, we compile photometry from the 
PS1 detections, our optical observations, and the circulars to 
construct a light curve for SN 2010ay. 

2.3. Optical Observations 

We obtained an optical spectrum (~3000-l 1000 A) of 
SN 2010ay on April 1 UT, from the ISIS blue arm instrument 
of the 4.2 m WHT at the Roque de Los Muchachos Observa- 
tory. The spectrum was taken at the parallactic angle, and the 
exposure time was 1800 s. We obtained a second 1800 s opti- 
cal spectrum (~3600-9600 A) using the Gemini Multi-Object 
Spectrograph (GMOS) on the 8.1 m Gemini North telescope on 
2010 April 11.4 UT. We employed standard two-dimensional 
long-slit image reduction and spectral extraction routines in 
IRAF. 21 We do not apply a correction for atmospheric differen- 
tial refraction, because the displacement should be <0?5 at the 
airmass of the observations, «H.0. 

In both our Gemini and WHT spectra, broad absorption 
features associated with the SN are clearly detected in addition 
to narrow emission lines typical of star-forming galaxies. 
We distinguish the host galaxy emission from the continuum 
dominated by the highly broadened SN emission by subtracting 
a high-order spline fit to the continuum. Both SN and host 


21 IRAF is distributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science Foundation. 
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Table 3 

Emission-line Fluxes Measured for the Host Galaxy of SN 201 Oay 


Emission Line 

Flux 

(10 -16 erg s^ 1 cm -2 ) 

[On] 443726, 3729 

52 ± 1 

H<5 

4.5 ± 0.2 

H y 

9.9 ±0.1 

[O hi] A4363 

0.7 ±0.1 

H ft 

24.5 ± 0.1 

[O hi] 44959 

30.4 ± 0.1 

10 hi] 45007 

90. 1 ± 0.2 

[Nn] 46548 

2.2 ±0.1 

H« 

84.5 ± 0.1 

[Nn] 46584 

6.33 ± 0.06 

[S 11 ] 46717 

7.46 ± 0.06 

[S 11 ] 46731 

5.59 ± 0.06 


Notes. All fluxes have been measured from our Gemini 
spectrum. No reddening correction has been applied. There is 
an additional systematic uncertainty in the flux measurements 
of -10% due to flux calibration. 

galaxy spectral components are shown in Figure 3. As illustrated 
in the figure, the broad, highly blended spectral features of 
SN 2010ay resemble those of the Type Ic-BL SN 2010bh 
(associated with GRB 100316D) at a similar epoch (Chornock 
et al. 2010). In particular, the broadening and blueshift of the 
feature near 6355 A are similar for SN 2010ay and SN 2010bh 
and are broader and more blueshifted than in SN 1998bw at a 
comparable epoch. We discuss the comparison between these 
two SNe further in Sections 3.3 and 6.2. 

Additionally, we obtained 60 s gri -band images of SN 2010ay 
on 2010 April 11.4 UT using GMOS. The data were reduced 
using the gemini package in IRAF, and photometry was 
performed using the standard GMOS zero points. 22 We measure 
that [g, r, i] = [18.90, 18.32, 18.04] ± 0.1 mag. 

Imaging photometry is not available at the epoch of our WHT 
spectrum. However, the spectrum was flux-calibrated against 
observations of the standard star Feige 34, which was observed 
the same night and at approximately the same low airmass as the 
SN. For the observations of both the standard star and SN, the 
slit was placed at the parallactic angle. We perform synthetic 
photometry on the spectrum to extract the flux at the central 
frequency of the R band (6527 A) and find R = 18.2 ± 0.2 mag. 
We then subtract the host galaxy flux numerically. 

2.3.1. Host Galaxy Features 

We measure fluxes of the narrow emission lines from the host 
galaxy in our Gemini spectrum, as reported in Table 3. We fit a 
Gaussian profile to each narrow emission line; for nearby lines 
such as [N 11] and Ha, we fit multiple Gaussians simultaneously. 
We model the local continuum with a linear fit to 20 A regions on 
either side of each line. We estimate uncertainties in quantities 
derived from the line fluxes by Monte Carlo propagation of the 
uncertainties in the flux measurement. 

The host galaxy is significantly reddened as evidenced by the 
flux ratio of Ha to H /3, «3.45 ± 0.02. We infer E(B - V) = 
0.2 mag (Ay = 0.6 mag), as measured from the Balmer 
decrement in our Gemini spectrum, assuming Ry = 3. 1, Case B 
recombination (Osterbrock & Ferland 2006), and the reddening 
law of Cardelli et al. (1989). This is similar to the value derived 


22 http://www.gemini.edu/?q=node/10445 
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4000 5000 6000 7000 8000 9000 

Rest wavelength (A) 


Figure 3. Optical spectra of SN 2010ay from Gemini (24 days after R-band peak) and the WHT (14 days). The spectrum is plotted decomposed into SN (above, 
with narrow lines clipped) and host galaxy (below, from Gemini, with spline-fit subtracted) components for clarity. The spectrum of SN 2010bh from Chomock et al. 
(2010) is given in black for comparison, at 21.2 days after the GRB 100316D trigger (~10.0 days after R-band peak; Cano et al. 201 la). The spectrum of SN 1998bw 
at +19 days from Patat et al. (2001) is also plotted. Both are transformed to the redshift of SN 2010ay. The SNe are shifted in flux for clarity. In the lower plot, relevant 
host galaxy emission lines are marked with a red line and labeled. 

(A color version of this figure is available in the online journal.) 


from the SDSS DR8 nuclear fiber spectroscopy line fluxes for 
the host galaxy ( E(B — V) = 0.2 mag). The value reported in 
Modjaz et al. (2010) was also similar: E(B — V) = 0.3 mag. 
Furthermore, we note that the color (B — V = 0.78 ± 0.07 mag) 
as reported by Prieto (2010) at 4 days after R-band peak is 
significantly redder than SN Ib/c color curve templates (Drout 
et al. 2011), further supporting a non-negligible host galaxy 
extinction. 

2.4. Radio Observations 

We observed SN 2010ay with the EVLA (Perley et al. 
2009) on three epochs, 2010 March 26, 2010 April 29, and 
2011 May 7. All EVLA observations were obtained with 
a bandwidth of 256 MHz centered at 4.9 GHz. We used 
calibrator J 1 22 1 +28 1 3 to monitor the phase and 3C 286 for flux 
calibration. Data were reduced using the standard packages of 
the Astronomical Image Processing System. We do not detect 
a radio counterpart to SN 2010ay in these observations and 
place upper limits of F v < 46, 42, 30 /iJy (3a) for each epoch, 
respectively, corresponding to upper limits on the spectral 
luminosity spanning L v < (3. 6-5.5) x 10 27 ergcm 2 s '. 

As shown in Figure 4, these limits are comparable to the 
peak luminosities observed for ordinary SNe Ib/c (Berger 
et al. 2003a; Soderberg 2007; Chevalier & Fransson 2006; 
Soderberg et al. 2010, and references within) and a factor of 
10 2 — 10 3 less luminous than the radio afterglows associated 
with GRBs 020903, 030329, and 031203 at early epochs 
(Berger et al. 2003b; Soderberg et al. 2004a, 2004b; Frail 


et al. 2005). In comparison with the radio luminosities observed 
for the relativistic SNe 1998bw (Kulkarni et al. 1998) and 
2009bb (Soderberg et al. 2010), SN 2010ay is a factor of 
>10 less luminous. The peak radio luminosity observed for 
GRB 100316D was a factor of a few higher than the first 
EVLA non-detection of SN 2010ay (E. Margutti et al. 2012, in 
preparation). The only relativistic explosion with detected radio 
emission below our EVLA limits is the weak and fast-fading 
XRF 060218 (Soderberg et al. 2006a). 


3. INITIAL CONSTRAINTS 
3.1. Light-curve Modeling 

We construct an R-band light curve for SN 2010ay using the 
observations described in Section 2.3 (Table 2). We convert the 
ip i - and r'-band data points to the R band using the fiducial 
light-curve method of Soderberg et al. (2006b), assuming the 
unextincted i ' — R and r' — R colors observed for the Type Ic-BL 
SN 1998bw at the appropriate epoch (Galama et al. 1998). 
Photometry for the unfiltered CRTS images was reported by 
Drake et al. (2010) after transformation to the synthetic V hand 
(A. J. Drake 201 1, private communication); we therefore convert 
to the R band assuming the V— R color of SN 1998bw at the 
appropriate epoch. For the Pan-STARRS 1 photometry, the host 
galaxy flux was subtracted using a template image. For all other 
photometry, we have subtracted the flux of the host galaxy 
numerically assuming the magnitude reported in Table 1 . A total 
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Figure 4. EVLA upper limits for SN 2010ay (black arrows) are compared with 
off-axis GRB afterglow light-curve models (black curves; standard parameters: 
Ek & 10 51 erg, A * = 1, € e = eg = 0.1, 0j = 5°, and viewing angles 
of 30°, 60°, 90°) and the observed radio light curves for ordinary SNe Ib/c 
(gray; Soderberg 2007, and references within) and the radio afterglows of all 
GRB-SNe within z ^ 0.25. SN 2010ay is a factor of 10 2 — 10 3 less luminous 
than XRF 020903 (orange; Soderberg et al. 2004a), GRB 030329 (blue; Berger 
et al. 2003b; Frail et al. 2005), and GRB 031203 (Soderberg et al. 2004b). 
Relativistic, engine-driven SNe 1998bw (red; Kulkarni et al. 1998) and 2009bb 
(dark blue; Soderberg et al. 2010) are a factor of 10 more luminous than the 
SN 2010ay limits on a comparable timescale, while XRF 060218 lies a factor 
of a few below the limits. We constrain the radio counterpart to be no more 
luminous than XRF 0602 1 8 and comparable to the peak luminosities of ordinary 
SNe Ib/c. 

(A color version of this figure is available in the online journal.) 

(Galactic + host) reddening of E(B — V) = 0.2 mag has been 
assumed (see Section 2.3.1). 

To estimate the explosion date of SN 2010ay, we have fit an 
expanding fireball model to the optical light curve (Figure 2), 
following Conley et al. (2006). In this model, the luminosity 
increases as 


We derive an explosion date to of 2010 February 21.3 ± 1.3. 
Here, we have assumed an index n = 2. This suggests that the 
PS1 37r z'pi - band pre-discovery detection image of SN 2010ay 
was taken ~4 days after the explosion. Observations by the 
PS1 survey have therefore provided a valuable data point for 
estimating the explosion date and also for constraining the 
rise time of SN 2010ay, as well as other nearby SNe (e.g., 
SN 201 lbm; Valenti et al. 2012). 

In Figure 2, we compare the light curve of SN 2010ay to 
the SN Ib/c light-curve template of Drout et al. (2011) after 
stretching by a factor of (1 + z). The template provides a 
reasonable fit to the optical evolution of SN 2010ay. Fitting the 
template to our photometry, we derive (reduced x 2 = 0.3) a date 
of R-band peak of 2010 March 18 ± 2 UT (MJD 55,273 ± 2) 
and an R-band peak magnitude of Mr ~ —19.7 mag before 
extinction correction. As discussed in Section 2.3.1, based on 
the Balmer decrement observed for the host galaxy emission 
lines, we assume an extinction of E(B — V) = 0.2 mag. 
Applying this correction, the peak absolute magnitude is Mr ~ 
—20.2 ± 0.2 mag. We note that this fitted value is ~0.2 mag 
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fainter than that estimated from the data point near peak. Here, 
the uncertainty is dominated by the template fitting. 

Regardless of the extinction correction, SN 2010ay is more 
luminous than all the 25 SNe Ib/c in the sample of Drout 
et al. (201 1), except for SN 2007D (Mr —20.65 mag, which 
was also significantly extincted: Ay ~ 1.0 mag). Assuming an 
intrinsic V— R color of zero at peak (e.g., 1998bw: Galamaetal. 
1998; Patat et al. 2001), SN 2010ay is also more luminous than 
any of the 22 GRB and XRF-producing SNe in the compilation 
of Cano et al. (2011b), all corrected for extinction, and is 
1.5 standard deviations from the mean luminosity. The peak 
magnitude is only ~ 1 mag below that of the Type Ic SN 2007bi 
(Mr = —21.3 ±0.1 mag), which Gal- Yam et al. (2009) report 
as a candidate pair-instability SN. 

3.2. Large Nickel Mass for SN 20 Way 

We use the available photometry for SN 2010ay discussed 
above to derive the mass of 56 Ni required to power the optical 
light curve under the assumption that the emission is powered 
by radioactive decay. Using the relation between M® and Mr 
found by Drout et al. (201 1), log(MNi) = (— 0.41Mr — 8.3) M Q , 
we estimate that SN 2010ay synthesized a nickel mass of 
M Ni = 0.9^0*! M q . We have estimated the uncertainty in the M« 
estimate by propagation of the uncertainty in the template fitted 
peak magnitude — systematic uncertainties are not included. If 
we instead adopt the most luminous individual data point in the 
light curve as the peak value, rather than the smaller peak value 
from template fitting, we find M Ni ~ 1 .2 M 0 . 

The Mn; estimate for SN 2010ay is larger than that of all 
but one (SN 2007D) of the 25 SNe Ib/c in the Drout et al. 
(201 1) compilation and significantly larger than the estimate for 
GRB-SN 2010bh, M Ni = 0.12 ± 0.01 M 0 (Cano et al. 2011a). 
On the other hand, Mni of SN 2010ay is at least 3 times smaller 
than for SN 2007bi (M Ni % 3.5-4.5M 0 ; Young et al. 2010). 
A pair-instability SN should produce M^ > 3 M 0 (Gal-Yam 
et al. 2009). 

3.3. Unusually High Velocity 

As illustrated in Figure 3, the broad, highly blended spectral 
features of SN 2010ay at the time of the WHT observations 
(14 days after R-band peak; see Section 3.1) resemble those of 
SN 2010bh at a similar time (10.0 days after peak; Chornock 
et al. 2010). In particular, the blueshift of the feature near 6355 A 
is larger than in SN 1998bw and more similar to SN 2010bh. 
This feature is commonly associated with Si ii A6355 (e.g., Patat 
et al. 2001). However, this feature has two clearly detectable 
absorption minima in SN 2010ay, but not in SN 2010bh. This 
could be due to increased blending in SN 2010bh or the absence 
of contaminating lines. The red ends of the SN 2010ay and 
SN 2010bh spectra (rest wavelength > 7500A) have similar 
P Cygni features, but the emission and absorption components 
in SN 2010bh are each blueshifted by ~200A more than in 
the spectrum of SN 2010ay. Chornock et al. (2010) attribute 
this feature to the Can near-IR (NIR) triplet, with a gf- 
weighted line centroid of 8479 A, and find a velocity that is 
high but consistent with the early-time velocity of Sin a 63 5 5 
((30-35) x 10 3 kins" 1 ). 

We measure the absorption velocity from the minimum of the 
Sin /. 63 5 5 absorption feature (t>si). We smooth the spectrum 
using an inverse-variance-weighted Gaussian filter (Blondin 
et al. 2006, with dX/X = 0.005) and measure the minimum 
position of the redmost component of the absorption profile. 
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Figure 5. Comparison of the time-evolving absorption velocity of SN 201 Oay 
and other SNe Ic-BL (red) and engine-driven explosions (blue) from the 
literature. For each SN, we fit a power law of the form t>si = «i P (f/30)“, 
where t is the time in days, i;^ 1 is the velocity at 30 days after explosion (dashed 
vertical line), and a is the velocity gradient. The velocities for SNe 1997ef, 
2003dh, and 20031w are from Mazzali et al. (2006), as determined by spectral 
modeling. The velocities for all other SNe are measured from the Si n A6355 
feature as follows: SN 2007ru is from Sahu et al. (2009); SNe 1998bw, 2006aj, 
and 2010bh are from Chomock et al. (2010), from spectra in references therein; 
SN 2007bg is from Young et al. (2010); SNe 2002ap, 2009bb, and 2003jd are 
fromPignata et al. (2011), and references therein; and SN 2010ay is from Prieto 
(2010) and this paper. See Table 4 for uncertainties. 

(A color version of this figure is available in the online journal.) 


Sanders et al. 


Table 4 

Velocity Evolution of SNe Ic-BL 


SN 

30 

u Si 

a 

SNe Ic-BL 

1997ef 

6 ± 2 

-0.8 ± 0.4 

2002ap 

5 ± 2 

-1.5 ±0.7 

2003jd 

10 ± 1 

-0.5 ± 0.2 

2007bg 

9 ± 2 

-0.3 ± 0.2 

2007m 

10.3 ± 0.7 

-0.5 ±0.1 

201 Oay 

21 ± 2 

-0.4 ± 0.3 

Engine-driven SNe Ic-BL 

1998bw 

10.3 ± 0.7 

-0.86 ± 0.08 

2003dh 

12.0 + 0.9 

-0.9 ±0.1 

20031w 

10+ 1 

-0.8 ± 0.2 

2006aj 

15 + 3 

-0.3 ± 0.2 

2009bb 

10.7 ± 0.4 

-0.84 ± 0.08 

201 Obh 

24 + 3 

-0.22 ± 0.07 


Notes. To the velocity measurements for each SN, we have fit a 
power law of the form usi = ttg| ) (f/30) Q ', where t is the time since 
explosion in days and is the velocity at 30 days in units of 
10 3 knis 1 . These power-law fits are illustrated in Figure 5. We adopt 
a conservative error of 10% on individual velocity measurements in 
order to compare the allowable range of fitted parameters for different 
SNe and note that this uncertainty is likely dominated by systematic 
effects. Uncertainties in the explosion date for SNe without detected 
GRBs vary due to availability of early-time photometry; we adopt 
the following conservative uncertainties: 7 days (1997ef; Hu et al. 
1997), 7 days (2002ap; Gal- Yam et al. 2002), 6 days (2003jd; Valenti 
et al. 2008), 14 days (2007bg), 3 days (2007ru; Sahu et al. 2009), 
1 day (2009bb; Pignata et al. 2011), and 2 days (2010ay; this paper). 


The blue component of the absorption profile shifts blueward 
over time, suggesting that it is produced by a combination of 
ions separated by several 10 3 km s _1 , such as Hei >,5876 and 
Nai D, whose relative optical depth changes with time. 

The absorption velocity inferred from the Sin 16355 feature 
is ~2 times faster than that of SN 1998bw at similar times, and 
more similar to that of SN 2010bh (Figure 3). For SN 1998bw, 
Patat et al. (2001) measured ~10 x 10 3 km s” 1 at +13 days. 
For SN 201 Obh, Chornock et al. (2010) measured usi ~ 
26 x 10 3 km s _ 1 at +10.0 days after explosion. Prieto (2010) 
reported a velocity of t>si ~ 22.6 x 10 3 km s -1 from the 
Si ii 16355 feature in a spectrum of SN 2010ay taken at +4 days; 
Prieto (2010) does not discuss the details of their methodology 
for the velocity measurement. From our [WHT, Gemini] spectra 
taken [+14, +24] days after /(-hand peak (see Section 3.1), we 
estimate i>si ~ [19.2, 18.3] x 10 3 km s _1 . 

In addition to the broadening of the spectral features and the 
blueshift of the Sin 16355 line, additional lines of evidence 
suggest a high photospheric expansion velocity for SN 201 Oay. 
We measure Vsi ~ [2 1 . 7 , 20 . 1 ] x 10 3 km s 1 from the Ca ii NIR 
triplet on the smoothed [WHT, Gemini] spectra, relative to a 
line center at 8479 A. This is within a few x 10 3 km s -1 of the 
u si we measure from Si ii 16355 at these epochs. Furthermore, 
we do not detect the broad emission “bump” near 4500 A in 
either of our spectra. This feature was also absent in SN 201 Obh 
but was identified in SN 2003dh, SN 2006aj, and several Ic-BLs 
not associated with GRBs and normal SNe Ic; Chornock et al. 
(2010) suggest that the absence of this feature indicates a high 
expansion velocity if it is due to blending of the iron lines to the 
blue and red of 4500 A. 

We compare the late-time expansion velocity of SN 201 Oay 
to other SNe Ic-BL and GRB-SNe with detailed, multi-epoch 


velocity measurements from the literature in Figure 5. We 
note that the velocities for SNe 1997ef, 2003dh, and 20031w 
plotted in the figure were estimated by Mazzali et al. (2006) 
via spectral modeling, rather than measured directly from the 
minimum of the Sin 16355 absorption feature; however, these 
velocities should be equivalent to within a few 10 3 km s -1 , 
as the minimum of the Si ii feature is typically well fit by the 
photospheric velocity of the spectral models (see, e.g., Mazzali 
et al. 2000; Kinugasa et al. 2002). In this figure, we also fit 
power-law gradients to the time evolution of the velocity of 
these SNe with the form usi = ^(//SO)®, where t is the time 
since explosion in days and is the velocity at 30 days in 
units of 10 3 km s . These parameters are listed in Table 4. 
Figure 5 illustrates that most SNe are well described by a single 
power law. 23 However, due to lack of late-time spectroscopy, the 
Ug? measurement amounts to an extrapolation for some objects 
(particularly 2006aj), and contamination from different ions or 
detached features will add uncertainty to velocities measured 
from the Sin >,6355 feature. 

SNe 201 Oay and 201 Obh share high characteristic velocities 
at 30 days after explosion and velocity gradients that are low 
relative to other broad-lined SNe Ic with and without associated 
GRBs. For SN 201 Oay, vf* = 21 ± 2 is 2-4 times larger than 
for other SNe Ic-BL without associated GRBs (u^ = 6 ± 2 
for 1997ef, 5 ± 2 for 2002ap, 10 ± 1 for 2003jd, 9 ± 2 for 
2007bg, and 10.3 ± 0.7 for 2007m) and is similar to the 
GRB-SN 201 Obh (u 3 ° = 24 ± 3). No other GRB-SN or SN Ic- 
BL has ul? > 15. The SNe Ic-BL and GRB-SNe with the 


23 A break at ~2 x 10 4 km s 1 appears to exist for SN 1998bw at ~16 days, 
as noted by Kinugasa et al. (2002). 
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most shallow velocity gradients among these 12 objects have 
a < -0.5; they are SN 2006aj (a = -0.3 ± 0.2), SN 2007bg 
(a = — 0.3±0.2), SN 2010ay (a = — 0.4±0.3), and SN 2010bh 
(i a = —0.22 ± 0.07). This places SN 2010ay in the company of 
two GRB-SNe in having a slowly evolving absorption velocity 
and SN 2007bg (whose unusually fast decline rate distinguishes 
it from other SNe Ic-BL; Young et al. 2010). The velocities 
of SNe 2010ay and 2010bh, respectively, decline about 2 and 4 
times more slowly than the other SNe Ic-BL (mean and standard 
deviation: a = —0.8 ± 0.4) and about 1.5 and 3 times more 
slowly than the other GRB-SNe ( a = —0.7 ± 0.2). The slow 
Sin k6355 absorption velocity evolution of SN 2010ay at late 
times resembles the slow evolution of the Fen lines of the 
spectroscopically normal Type 1c SNe 2007gr and 201 lbm at 
late times (Valenti et al. 2012). 

Given the high peak luminosity of SN 2010ay (Section 3.1), 
we also consider the velocity of the candidate pair-instability 
SN 2007bi. Velocity measurements for SN 2007bi are only 
available at late times (>50 days; Young et al. 2010). Fitting 
to these late-time Sin 16355 velocity measurements, we find 
that SN 2007bi has a characteristic velocity —3 times smaller 
than 2010ay (u|° = 8 ± 2) and the late-time velocity gradient is 
—2 times more shallow ( a = —0.2 ± 0.2). 

3.4. Ejecta Mass and Energy 

We use the scaling relations provided by Drout et al. (2011), 
based on the original formalism of Arnett (1982) and modified 
by Valenti et al. (2008), to derive the total mass of the ejecta and 
the kinetic energy 


Me j = 0.81 ^ 


USi 


8 d / V 10, 000 km s- 1 


Mq 


( 2 ) 


E k — 0.5 — 


VSi 


8 d / \ 10, 000 km s- 


— 1 x 10 51 erg. (3) 


We assume the fitted peak magnitude for SN 2010ay (M s = 
—20.2 ± 0.2; Section 3.1), the absorption velocity we measure 
from our WHT spectrum at 14 days after maximum light 
(i>si = 19.2 x 10 3 km s' 1 ; Section 3.3), and a characteristic time 
(light-curve width) consistent with the data and the mean value 
from the Drout etal. (2011) sample of SNe Ic-BL (r c = 14 days). 

Using these values, the total mass ejected was M e j ~ 
4.7 M q , and the total kinetic energy of the explosion was 
Ek — 10.8 x 10 51 erg. Hereafter we refer to the definition 
£*,51 = £*/10 51 erg. 

The systematic uncertainties associated with this modeling 
dominate the statistical uncertainties. In particular, the models 
rely on the assumptions of homologous expansion, spherical 
symmetry, all 56 Ni centralized at the center of the ejecta, 
optically thick ejecta, and constant opacity. 

We note that an earlier measurement of the absorption velocity 
is preferred for optical modeling. Since we have argued that 
SN 2010ay and SN 2010bh have similar characteristic velocities 
( v s?), if we instead adopt a higher velocity of 25,000 km s” 1 as 
measured for SN 2010bh by Cano et al. (2011a), we estimate 
Me j « 6.1 Mq and Ek . 51 — 23.9 for SN 2010ay. 

The M e j and E K 51 of SN 2010ay are consistent with the mean 
for SNe Ic-BL in the Drout et al. (2011) sample (4.7* 2 f 3 8 Mq 
and 1 U4, respectively), because the authors assumed a velocity 
(usi = 2 x 10 4 km s” 1 ) similar to the late-time velocity we 
measure. For comparison, SN 2010bh had a total ejecta mass of 


—2 Mq and a total kinetic energy of E K 5 j — 13 (Cano et al. 
2011a). 

The ratio of Ni to total ejecta mass is —0.2 for SN 2010ay, 
significantly higher than the values typical of SNe Ic-BL 
and GRB-SNe. For comparison, the ratio is just —0.05 for 
SN 2010bh. Adopting the values derived from bolometric light- 
curve modeling by Cano et al. (2011a), the M® and MNi/M e j 
ratios for other GRB-SNe are —0.5 M 0 and -0.06—0.22 
(1998bw), ~0.4 M 0 and —0.08 (2003dh), — 0.15M o and 
-0.07-0.1 (2006aj), and ~0.2 M Q and -0.06 (2009bb). This 
ratio for SN 2010ay is larger than that of all but 4 of the 25 SNe 
ofDroutetal. (201 1): the Type Ic SNe 2004ge (MNi/M e j — 0.4) 
and 2005eo (M^/ M e j — 0.2), the Type lb SN 2005hg (—0.4), 
and the Type Ic-BL SN 2007D (-0.6). 

The large value of Mn; we estimate for SN 2010ay raises 
the question of whether a process other than Ni decay may be 
powering its light curve (e.g., Chatzopoulos et al. 2009). An 
independent test of the physical process powering the light 
curve is the decay rate of the late-time light curve, which 
should be «0.01 mag day” 1 for SNe powered by radioactive 
decay of 56 Co. For SN 2007bi, Gal- Yam et al. (2009) derive 
Mm = 3.5 M 0 from the measured peak magnitude and find that 
the late-time light curve is consistent with the decay rate of 56 Co. 
While the Pan-STARRS 1 3 n survey also observed the field in 
2011 March, the SN was not detected in our subtracted images, 
and the limits are not constraining in the context of 56 Co decay 
(see Section 2.2 and Table 2). Another possible process is a 
radiation-dominated shock that emerges due to interaction with 
an opaque circumstellar medium (CSM), as has recently been 
proposed by Chevalier & Irwin (2011) for the ultraluminous 
SNe 2006gy and 2010gx (Pastorello et al. 2010; Quimby et al. 
2011). However, while this class of ultraluminous objects shares 
some spectroscopic similarities to SNe 1c (Pastorello et al. 
2010), they show peak luminosities —4-100 times higher than 
SN 2010ay (Chomiuk et al. 2011; Quimby et al. 201 1). 

4. CONSTRAINTS ON RELATIVISTIC EJECTA 

We use our EVLA upper limits for SN 2010ay spanning 
At ~ 29-433 days to constrain the properties of the shock wave 
and those of the local circumstellar environment. The radio 
emission from SNe Ib/c and GRBs is produced by the dynami- 
cal interaction of the fastest ejecta with the surrounding material 
(Chevalier 1982). The kinetic energy of the ejecta is converted, 
in part, to internal energy of the shocked material, which itself 
is partitioned between relativistic electrons (e c ) and amplified 
magnetic fields (e ft ). Following the breakout of the shock wave 
from the stellar surface, electrons in the environment of the 
explosion are shock-accelerated to relativistic velocities with 
Lorentz factor y e and distributed in a power-law distribution 
characterized by N(y e ) oc y e p . Here, p characterizes the elec- 
tron energy index. The particles gyrate in amplified magnetic 
fields and give rise to non-thermal synchrotron emission that 
peaks in the radio and mm-bands in the days to weeks fol- 
lowing explosion with observed spectral index F v oc u” (p ” | 7 2 . 
At lower frequencies the emission is suppressed due to syn- 
chrotron self-absorption, which defines a spectral peak, v p 
(Chevalier 1998). 

The dynamics of the shock wave determine the evolution 
of the synchrotron spectrum and, in turn, the properties of the 
observed radio light curves. In the case of SNe Ib/c, there are 
three primary scenarios for the dynamical regime of the ejecta 
depending on the shock velocity, v = Pc (associated Lorentz 
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factor, T): (1 ) non-relativistic ( v 0.2c) free expansion as in the 

case of ordinary SNe Ib/c (Chevalier 1998); (2) a decoupled and 
relativistic (r ~ 10) shell of ejecta that evolves according to the 
Blandford-McKee solution for several months (Sari et al. 1998) 
before transitioning to the Sedov-Taylor regime (Frail et al. 

2000; this is the standard scenario for typical GRBs); and (3) a 
sub-energetic GRB with trans-relativistic velocity (/3T < 3) that 
bridges the free-expansion and Blandford-McKee dynamical ^ 
regimes (e.g., SN 1998bw; Kulkarni et al. 1998; Li & Chevalier % 
1999). ~ 

We consider our EVLA upper limits in the context of these 
three models below. For shock velocities of v > 0.2c, e e & 0. 1 is c 
reasonable (Soderberg et al. 2005; Chevalier & Fransson 2006). 

We further assume equipartition, e e = ( b = 0.1. We adopt 
a free-expansion model for both the non-relativistic ordinary 
SN Ib/c case and the sub-energetic, trans-relativistic GRB 
scenario. As shown by Li & Chevalier (1999), a free-expansion 
model is still reasonable in the trans-relativistic regime (case 3, 
see above). 

4.1. Freely Expanding Shock Wave 

In the free-expansion scenario, a shock discontinuity sepa- 
rates the forward and reverse shocks, located at the outer edge of 
the stellar envelope. The bulk ejecta is in free expansion, while 
the thin layer of post-shock material is slightly decelerated, 

R oc r 0 9 (Chevalier & Fransson 2006). At a given frequency, 
the bell-shaped light curves of the SN synchrotron emission may 
be described as (Chevalier 1998) 



Figure 6. Region of energy-velocity space ruled out (red) by our EVLA 
observations for on-axis ejecta under the assumption of a free-expansion model. 
(A color version of this figure is available in the online journal.) 

properties of M = 1 0 -5 M Q yr - 1 and a progenitor wind velocity 
of v w = 10 3 kms -1 : 


L v 


1.582 x L v , p 



[1 _ e -(A 


(4) 


where L v p is the flux density at the spectral peak at epoch 
t p . Assuming an electron index of p ~ 3, consistent with 
radio spectra of SNe Ib/c (Chevalier & Fransson 2006), the 
exponents are a & 2.3 and h ~ 1.3. The time-averaged shock 
wave velocity is v ^ R/At, where R is the shock wave radius 
defined as 


R ^ 2.9 x 10 16 | — 

, e B 


- 1/19 


10 28 erg s -1 Hz~ 


9/19 


V5GHz/ 


cm. 


(5) 


Here, we make the assumption that the radio-emitting region 
is half of the total volume enclosed by a spherical blast wave. 
Next, we estimate the internal energy, E, of the radio-emitting 
material from the post-shock magnetic energy density, E ss 
B 2 R 3 /12cb, where we maintain the assumption of equipartition. 
As shown by Chevalier (1998), the amplified magnetic field at 
peak luminosity is also directly determined from the observed 
radio properties, 


B 




J V,p 


10 28 erg s -1 Hz 


- 2/19 


(6) 


Finally, the mass-loss rate of the progenitor star, M, may be 
derived from the number density of emitting electrons. Here, 
we normalize the wind profile according to p oc Ar ~ 2 and 
A* = A/5 x 10 11 g cm - . This normalization of A* implies that 
an A* of 1 corresponds to typical Wolf-Rayet progenitor wind 


0.15 


0A 

v p 

5 GHz 


- 8/19 


A t V 


10 28 erg s -1 Hz 


-l 


- 4/19 


10 days. 


(7) 


where we assume a shock compression factor of ~4 and a 
nucleon-to-electron ratio of 2. 

We built a two-dimensional grid of fiducial radio light curves 
according to Equation (4) in which we vary the parameters L „ 
and v p over a reasonable range of parameter space, bounded by 
t p ~ [1,3000] days and F v p ~ 0.04-1000 mJy. We identify the 
fiducial light curves associated with a radio luminosity higher 
than the EVLA upper limits for SN 2010ay at each epoch as 
these are excluded by our observations. We extract the physical 
parameters associated with these excluded light curves (R, B, 
E, A*) to define the parameter space excluded by our radio 
observations. The parameter space for v p and F v p is bounded 
by the respective values for which the model exceeds relativistic 
velocities, ySF ~ few. 

As shown in Figure 6, our deep EVLA limits enable us 
to rule out a scenario in which there is copious energy, 
E > 10 48 erg, coupled to a relativistic outflow, in this two- 
dimensional E — v parameter space. The excluded region in- 
cludes GRB-SNe 1998bw and 060218, as well as the relativistic 
SN 2009bb. It does not exclude the standard scenario in which 
a small percent of the energy is coupled to fast-moving ma- 
terial within the homologous outflow, as is typically observed 
for ordinary SNe Ib/c (E ~ 10 47 and v ~ 0.2c; Soderberg 
et al. 2010). 

Next we consider the effects of circumstellar density since 
lower mass-loss rates produce fainter radio counterparts. As 
shown in Figure 7, the EVLA limits for SN 2010ay exclude 
the region of parameter space populated by SNe 1998bw and 
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Figure 7. Region of energy-mass-loss space ruled out (red) by our EVLA 
observations for on-axis ejecta under the assumption of a free-expansion model. 
(A color version of this figure is available in the online journal.) 

2009bb with mass-loss rates of A* ~ 0.1; however, the low- 
density environment of GRB 0602 1 8 lies outside of our excluded 
region due to its lower CSM density, A* ~ 0.01, which gives 
rise to a lower luminosity radio counterpart. 

4.2. Relativistic Ejecta 

In the case of relativistic deceleration the ejecta are confined 
to a thin jet and are physically separated from the homologous 
SN component. Deceleration of the jet occurs on a timescale 
of At R< (ZS51/A*) years in a wind-stratified medium (Waxman 
2004). On this same timescale, any ejecta components that were 
originally off-axis spread sideways into the observer’s line of 
sight. While the early EVLA limits constrain the properties 
of the on-axis ejecta according to the free-expansion model 
described above, the late-time EVLA upper limits constrain any 
radio emission from a GRB jet originally pointed away from 
our line of sight. 

For this scenario, we adopt the semianalytic model of 
Soderberg et al. (2006c) for off-axis GRB jets. In Figure 4, 
we compare the radio upper limits for SN 2010ay with the pre- 
dictions for an off-axis jet with standard parameters ( E K ~ 
10 51 erg. A* = 1, € e = €b = 0.1, and Oj = 5°) and assuming a 
viewing angle of 0 oa = 30°, 60° , or 90° . Our EVLA upper limits 
rule out all three of these model light curves. We derive the two- 
dimensional parameter space (energy and CSM density) that is 
excluded based on our EVLA upper limit at At ~ 1.2 years. 
We note that this model accommodates the full transition from 
relativistic to non-relativistic evolution. We built a collection of 
model light curves spanning parameter range A* Rs [0.01-100] 
and E ~ [10 49 -10 52 ] erg, maintaining the assumption of p = 3 
and e e = €b = 0.1. Here, we adopted a jet opening angle of 
9j = 5° and an off-axis viewing angle of 0 oa = 90° (the most 
conservative scenario). We note that van Eerten & MacFadyen 
(201 1 ) have developed off-axis GRB afterglow light-curve mod- 
els based on hydrodynamic simulations that reproduce the semi- 
analytic models presented here to within factors of a few (see 
also van Eerten et al. 2010). 



Mass Loss Rate (A.) 

Figure 8. Our EVLA observations at At 1 .2 years after explosion constrains 
the properties of a possibly associated off-axis GRB jet. Using our semianalytic 
model as described in Section 4.2, we assume partition fractions of = eg = 
0.1, Oj = 5°, p = 2.5, and an off-axis viewing angle of # oa = 90°. We are able 
to exclude the region of F.k — A . parameter space (red) associated with typical 
GRBs, i.e., Eg = 10 51 erg and A* = 1 (dashed black lines). 

(A color version of this figure is available in the online journal.) 

As shown in Figure 8, we are able to exclude the parameter 
space associated with typical GRBs, i.e., £ R< 10 51 erg (beaming 
corrected) and A* R* 1. GRBs with lower energies and densities 
are better constrained using the trans-relativistic formalism 
above. In conclusion, our radio follow-up of SN 2010ay reveals 
no evidence for a relativistic outflow similar to those observed 
in conjunction with the nearest GRB-SNe: however, a weak 
afterglow like that seen from XRF 060218 cannot be excluded. 

5. CONSTRAINTS ON AN ASSOCIATED GRB 

Given the estimate of the explosion date we have derived 
(Section 3.1), we have searched for gamma-ray emission that 
may have been detected by satellites. No GRBs consistent with 
SN 2010ay were reported in the circulars of the Gamma-ray 
Coordinates Network Circulars, but it is possible that bursts 
were detected below the instrument triggering thresholds. 

We next consulted the sub-threshold bursts from the Swift 
Burst Alert Telescope (BAT; Gehrels et al. 2004; Barthelmy 
et al. 2005) detected within the 6 days surrounding the explosion 
date estimate. We find that no gamma-ray emission was detected 
within 0°5 of the position of SN 2010ay by the BAT during this 
period. Given the sensitivity of the BAT, this corresponds to 
an upper limit on the gamma-ray flux of ~10~ 8 erg s -1 cm ~ 2 
(15-150 keV). However, the field of SN 2010ay was in the field 
of view of the instrument for only 106 ks during these 6 days, 
or ~20% of the duration. 

For complete temporal coverage, we have searched the 
records of the interplanetary network (IPN), which is sensitive to 
bursts with fluences down to ~6 x 10~ 7 erg cm -2 (25-150 keV) 
(50% efficiency limit; Hurley et al. 2010) and observes the 
entire sky with a temporal duty cycle close to 100%. An 
undetected, sub-threshold burst should have a fluence below 
this limit. Between 2010 February 21 and 25, inclusive, a total 
of 12 bursts were detected by the spacecraft of the IPN (Mars 
Odyssey, Konus-Wmr/, RHESSI, INTEGRAL (SPI-ACS), Swift- 
BAT, Suzaku, AGILE, MESSENGER, and Fermi (GBM)). Ten 
of them are confirmed bursts; they were observed by more than 
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one instrument on one or more spacecraft and could be localized. 
Two of them are unconfirmed bursts; they were observed by one 
instrument on one spacecraft ( Suzaku ). The total area of the 
localizations of the 10 confirmed bursts is ~0.58 x 4;r sr. This 
implies that about 0.58 bursts can be expected to have positions 
that are consistent with any given point on the sky simply by 
chance (i.e., within the 3a error region), and indeed none of 
the bursts in this sample have a position consistent with the SN 
position. 

These non-detections imply upper limits to the gamma-ray 
energy ( E y ) of a burst that may have been associated with 
SN 2010ay. The IPN non-detection indicates E y < 6 x 10 48 erg 
(25-150 keV), while the BAT non-detection indicates that the 
peak energy of the burst was <1 x 10 47 erg s -1 (15-150 keV) if 
the burst occurred while in the field of view of the instrument. 

We consider whether or not a hypothetical GRB associated 
with SN 20 1 Oay would have been detected by BAT or the IPN if it 
had characteristics similar to well-studied SN-associated GRBs. 
The isotropic prompt energy release of long GRBs is typically 
Ey t i S o ~ 10 52 erg; however, the prompt emission of the sub- 
energetic class of GRB-SNe can be several orders of magnitude 
fainter (Soderberg et al. 2006a). GRB 980425/SN 1998bw 
had a peak gamma-ray luminosity of ~5 x 10 46 erg s' 1 
(24-1820 keV; Galama et al. 1998), which is a factor of 
two below our BAT limit, and E y y so ~ 5 x 10 47 erg (Pian 
et al. 2000), more than a factor of five below our limit. 
Neither satellite should have detected such a burst. In contrast, 
GRB 031203/SN 20031w had a peak gamma-ray luminosity of 
~lx 10 49 erg s -1 (20-200 keV) and a total isotropic equivalent 
energy of E y< j s0 = (4 ± 1) x 10 49 erg (Sazonov et al. 2004), 
about two orders of magnitude above the sensitivity of the BAT 
and twice the threshold of the IPN, respectively. GRB 030329/ 
2003dh was even more luminous, with E y i s0 ~ 7 x 10 49 erg 
(Hjorth et al. 2003). A burst like GRB 031203 or 030329 
should certainly have been detected by IPN, or the BAT if it 
occurred while the field of SN 201 Oay was in the field of view 
of the instrument. XRF 060218/SN 2006aj, an extremely long 
duration (At ~ 2000 s) event, had a peak luminosity observed by 
BAT of — 1 0 s erg cuT 2 s -1 (15-150 keV), corresponding to 
a peak emission of » 2.4 x 10 46 erg s” 1 given the redshift 
of the burst (z = 0.033), and a total isotropic equivalent 
energy of E Yt ; so = (6.2 ± 0.3) x 10 49 erg (Campana et al. 
2006). If such a burst was associated with SN 2010ay, its 
peak emission may have been a factor of four below the BAT 
sensitivity limit, although its total isotropic energy emission is 
an order of magnitude larger than our IPN limit for SN 201 Oay. 
Finally, the event whose host galaxy and SN properties seem 
most similar to SN 2010ay, GRB 100316D/SN 2010bh, had 
£y,iso > (5.9 ± 0.5) x 10 49 erg (Starling et al. 2011)— a full 
order of magnitude above our IPN limit. 

Another possibility is that prompt emission associated with 
SN 2010ay may have been too soft to be detected by the BAT 
or IPN. For example, the spectrum of XRF 060218 rose to a 
peak at 0.3-10 keV at ~985 s after triggering and then softened 
significantly thereafter. Even though the total emission of this 
burst is well above our IPN limit, it may have escaped detection 
if it was similarly soft. 

6. SUB-SOLAR HOST ENVIRONMENT METALLICITY 

We estimate the oxygen abundance of the host environment 
of SN 201 Oay from the strong nebular emission line fluxes 
measured from our Gemini spectrum (Table 3). At the distance 
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of the host galaxy, the 1" Gemini slit width corresponds to 
a physical size of 1.3 kpc. The properties we infer for the 
explosion site of SN 201 Oay represent a luminosity-weighted 
average over this radius. 

We employ several different oxygen abundance diagnos- 
tics in order to determine the metallicity of the host galaxy 
from its optical spectrum (e.g., Modjaz et al. 2011). From the 
03N2 diagnostic of Pettini & Pagel (2004, PP04), we derive 
a metallicity of log(0/H) = 8.19, or Z ~ 0.3 Z Q , adopt- 
ing the solar metallicity log(O/H) 0 + 12 = 8.69 from As- 
plund et al. (2005). Using the N2 diagnostic of PP04, we 
find log(0/H) + 12 = 8.26. Using the abundance diagnos- 
tic, R 2 3 = log([On] 73727 + [O in] 774959, 5007)/H/J, we 
find log(0/H) + 12 z94 = 8.49 (Zaritsky et al. 1994) and 
log(0/H) + 12 kdo 2 = 8.50 (Kewley & Dopita 2002). However, 
these /? 23 -based estimates are more sensitive to flux calibration 
and reddening correction. Moreover, there is a well-known bi- 
valued relationship between R 23 and oxygen abundance. The 
value R 23 = 0.873 ± 0.006 measured at the explosion site 
places it near the turnover point, but we assume that it lies 
on the upper branch based on its [Nil] /[On] ratio, following 
Kewley & Ellison (2008). The metallicity values we derive using 
the PP04 and KD02 diagnostics are approximately equiva- 
lent given the offset that exists between these two diagnostics 
(Kewley & Ellison 2008). These measurements are similar to the 
values reported by Modjaz et al. (2010) (log(0/H)+12 [PP04, 
KD02] = [8.2, 8.4]) for the SN 201 Oay host galaxy. The sta- 
tistical errors in our strong-line metallicity estimates are small 
(<0.01 dex), as determined by propagating the errors in the line 
flux measurement through the abundance calculation. However, 
for example, the representative systematic error for the PP04 
03N2 abundance diagnostic is ~0.07 dex, as determined by 
Kewley & Ellison (2008) via comparison to other strong-line 
abundance indicators. 

Fortunately, our detection of the weak [O in] 74363 auroral 
line (S/N ~ 6; Figure 3) allows us to derive an oxygen abun- 
dance via the “direct,” T e method. We employ a methodology 
similar to that used by, for example, Levesque et al. (2010a). We 
first derive the electron temperature (T e = (1 ,09±0.06) x 10 4 K) 
and density ( n e — 80 ± 20 cm 3 ) from the [O hi] and [S 11 ] line 
ratios using the temden task of the IRAF package nebular 
(Shaw & Dufour 1994), derive the 0 + temperature using the 
calibration of Garnett (1992), and finally estimate the 0 + and 
0 ++ abundances following Shi et al. (2006). The direct abun- 
dance, log(0/H) + 12 = 8.24 ± 0.08, is in good agreement 
with the PP04 03N2 value. The stated uncertainty reflects the 
propagation of the uncertainties for the line flux measurements. 
Indeed, the offset between these two diagnostics should be very 
small at this metallicity (Kewley & Ellison 2008). 

We estimate the star formation rate (SFR) of the host galaxy 
using the Ha relation of Kennicutt (1998). After correcting for 
host galaxy extinction, we measure the Ha luminosity from our 
Gemini spectrum (Table 3) and estimate SFR =1.1 M 0 yr _ 1 . 

6.1. Blue Compact Galaxy Host 

We compare the host galaxy of SN 201 Oay to the nearby 
galaxy population of the SDSS spectroscopic survey. The phys- 
ical properties of the host galaxy, SDSS J123527. 19+270402.7, 
are estimated in the MPA/JHU catalog. 24 The total (photomet- 
ric) galaxy stellar mass (M*) is given as 3. 6^' 3 x 10 s M 0 , 

24 http://www.mpa-garching.mpg.de/SDSS (described in Kauffmann et al. 
2003; Tremonti et al. 2004; Brinchmann et al. 2004; Salim et al. 2007, and 
updated for SDSS DR7). 
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the aperture-corrected SFR is 1 .0 +, (p Af 0 yr 1 , and the nuclear 
(fiber) oxygen abundance (0/H fl ) is log(0/H) + 12 = 8.58 /q q 2 3 
on the scale of Tremonti et al. (2004, T04). The specific star for- 
mation rate (SSFR) of the host galaxy is then ^2. 8^4 Gyr 1 . 
For consistency, we consider these values of Af*, the oxygen 
abundance, and the SFR for the host galaxy of SN 2010ay when 
comparing to other galaxies in the MPA/JFIU catalog. 

The oxygen abundance and SFR of the host galaxy of 
SN 2010ay listed in the MPA/JHU catalog are consistent with 
the values we derive in this paper (see also Kelly & Kirshner 
2011). The MPA/JHU catalog lists metallicities on the T04 
scale. Using the Kewley & Ellison (2008) conversion to the PP04 
scale, the T04 metallicity estimate corresponds to a metallicity 
of log(0/H) + 12 = 8.38, which is ~0.2 dex higher than the 
one we measure (log(0/H) +12= 8.19). However, there is a 
large (~0.2 dex) rms scatter between the PP04 03N2 and T04 
diagnostics at the regime of log(0/H)+ 12ppo4 ~ 8.2 (Kewley & 
Ellison 2008). The SFR in the MPA/JHU catalog is also in good 
agreement with the value we estimate from the Ha luminosity. 
Although our estimate does not include an aperture correction, 
the size of the Gemini slit (1) should encompass most of the star 
formation in the galaxy (Petrosian r = IT 355; Table 1). 

The mass-to-light ratio of the host galaxy of SN 2010ay is low 
compared to typical star-forming galaxies. To compare the host 
galaxy to the general galaxy population, we select a subset of 
the MPA/JHU catalog by requiring that estimates of Af*, SFR, 
and 0/H c be available, and we remove active galactic nuclei 
(AGNs) according to Kauffmann et al. (2003). We consider 
167,837 starbursting galaxies following these constraints. The 
host galaxy ranks in the [4th, 38th, 11th] percentile in [Af*, 
SFR, 0/H 0 ] among these galaxies. Among the selected galaxies 
with a stellar mass as low as the host galaxy, 23 the median 
and standard deviation of the Z?-band 26 absolute magnitude are 
— 15.8 ± 1.3 mag. With M B = —18.35 ± 0.05 mag, the host of 
SN 201 Oay is more luminous than other galaxies with a similar 
mass at the 2a level. The discrepancy is due to the blue color of 
the SN 2010ay host galaxy, which indicates a stellar population 
that is very young and therefore has a low stellar-mass-to-light 
ratio. Among the 1184 galaxies in the MPA/JHU catalog that 
meet the constraints above and have a color similar to the host of 
SN 2010ay (0.47 < u—r < 0.67, from SDSS fiber magnitudes), 
the host galaxy has typical properties, with [Af*, SFR, 0/H 0 ] in 
the [46th, 49th, 54th] percentile. 

Based on these properties, we classify the host galaxy of 
SN 2010ay as a luminous blue compact galaxy (BCG). BCGs 
span a large range in luminosity (—21 < M B < —12, where 
luminous BCGs have M B < —17) but are distinguished by 
their blue colors (B — V < 0.45), high SFR (1 < SFR < 
20 M q yr -1 ), and low metallicity (Z o /50 < Z < Z Q /2; Kunth 
& Ostlin 2000; Kong & Cheng 2002). The host galaxy of 
SN 2010ay has a luminosity ( M B = —18.35 ± 0.05), color 
(B — V = 0.11 ± 0.07), SFR (1.0!° 0 3 2 Af 0 yr -1 ), and metallicity 
(Z ~ 0.3 Z Q ) consistent with all these ranges. 27 


25 This subset is selected such that the host galaxy of SN 2010ay has the 
median mass: 1.39 < M* < 5.62 x 10 s Mq, /V s j m = 6, 978. 

26 We obtain B-band magnitudes by converting the Jr-corrected gri 
magnitudes given in the MPA/JHU catalog to BVR magnitudes using the 
transformation of Blanton & Roweis (2007). 

27 We note that a large fraction of luminous BCGs show evidence for 
disturbed morphologies or interaction with close companions (Garland et al. 
2004; Lopez-Sanchez et al. 2006), but we do not see evidence for a companion 
at the depth of SDSS images of the host galaxy of SN 2010ay. 
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6.2. Comparison to SN Ic-BL and GRB-SN Host Galaxies 

Our measurement of the metallicity from the Gemini spec- 
trum indicates that the explosion site of SN 2010ay is 
~0.5(0.2) dex lower in metallicity than the median SNe Ic 
(Ic-BL) in the sample of Modjaz et al. (2011). In that sample, the 
median PP04 03N2 metallicity measured at the explosion site 
of SNe Ic is log(0/H)+ 12 ^ 8.7 and for Ic-BL is ^8.4 dex, for 
12 and 13 objects, respectively. If instead the KD02 metallicity 
is used, the median of the sample is «8.9 dex for SNe Ic (13 
objects) and «8.7 dex for Ic-BL (15 objects), so the abundance 
of the SN 2010ay host galaxy is similarly low compared to the 
median. 

The metallicity of the environment of SN 2010ay is more 
similar to previously studied nearby GRB-SN progenitors. 
A metallicity identical to our measurement was measured 
at the explosion site of SN 2010bh (Levesque et al. 2011): 
log(0/H) +12 = 8.2. In the survey of Levesque et al. (2010a), 
and adding the measurement for SN 2010bh, the GRB-SN 
host galaxies have an average and standard deviation PP04 
03N2 metallicity of log(0/H) + 12 = 8.1 ± 0.1 on the PP04 
scale, which is consistent with the SN 2010ay environment. 
Among the 17 LGRB host galaxies surveyed in Savaglio et al. 
(2009), the average metallicity is somewhat lower, 1 /6 Z Q or 
log(0/H) + 12 ~ 7.9, but these are at an average redshift of 
z ~ 0.5 that is much higher than SN 2010ay. 

This evidence suggests that the host galaxy of SN 2010ay has 
chemical properties more consistent with LGRBs/GRB-SNe 
than SNe Ic-BL without associated GRBs; however, selection 
effects may mitigate this discrepancy. SNe found in targeted 
surveys of luminous galaxies have host galaxy properties biased 
toward higher metallicities, due to the luminosity-metallicity 
(L — Z) relation (Tremonti et al. 2004). LGRBs are found in 
an untargeted manner through their gamma-ray emission and 
therefore are not biased by this relation. 

SN 2010ay joins a growing list of SNe Ic-BL that have been 
discovered in low-metallicity host galaxies. Given the system- 
atic uncertainty in strong-line oxygen abundance diagnostics 
(~0.07 dex), we will consider host galaxies with metallicity 
Iog(0/H)ppo4 + 12 < 8.3 (Z < 0.4 Z Q ) to be in the low- 
metallicity regime of SN 2010ay. Among the 15 SNe Ic-BL 
(9 discovered by untargeted searches) in the surveys of Modjaz 
et al. (2008) and Modjaz et al. (2011), 4 were found in 
low-metallicity environments: SN [2007eb, 2007qw, 2005kr, 
2006nx] at log(0/H)ppo 4 + 12 = [8.26, 8.19, 8.24, 8.24], All 
of these SNe were discovered by untargeted searches. Young 
et al. (2010) measure the metallicity of the host galaxy of the 
broad-lined Ic SN 2007bg to be log(O/H)p P0 4 + 12 = 8.18, 
although this SN has light curve and spectral properties that dis- 
tinguish it from normal SNe Ic-BL (Section 3.3). Furthermore, 
Arcavi et al. (2010) find that SNe Ic-BL are more common in 
dwarf (Af,. ^ — 18) host galaxies, which the authors attribute to 
a preference for lower metallicities. 

The star formation properties of the host galaxy of SN 20 1 Oay 
also resemble the host galaxies of LGRBs. If we consider those 
galaxies in the MPA/JHU catalog with masses similar to the 
host galaxy of SN 201 Oay (as defined above), then the me- 
dian SFR and 0/H 0 of these galaxies is 0.13 Af 0 yr -1 and 
log(0/H 0 )+12 = 8.36, respectively. The host galaxy of 
SN 201 Oay is in the [96th, 77th] percentile for [SFR, O/HJ 
among these galaxies. This indicates that, while the host galaxy 
of SN 20 1 Oay falls within 1 a of the mass-metallicity (Af — Z) re- 
lation for star-forming galaxies, its SFR is extreme for its mass. 


12 


The Astrophysical Journal, 756:184 (16pp), 2012 September 10 



Host M b 

Figure 9. Plot of host galaxy metallicity vs. absolute B magnitude for 
SNe Ic-BL (red) and engine-driven explosions (blue). The L — Z relation of 
nearby star-forming galaxies is plotted as a solid line, with the 15th and 85th 
percentile boundaries of the galaxy distribution (dashed lines). Here, we have 
transformed the ^-corrected gri magnitudes from the MPA/JHU catalog to 
B band (Blanton & Roweis 2007) and converted the T04 metallicity values 
to the PP04 scale (Kewley & Ellison 2008), for the purpose of comparing it 
to metallicity measurements for SN host galaxies in the literature. The dot- 
dashed horizontal line is the divider between GRB-SNe and SNe Ic-BL host 
environments suggested by Modjaz et al. (2008). The host galaxy properties 
of GRB/SNe other than 2010ay are from the following references: Starling 
et al. 2011; Cano et al. 2011a (2010bh), Modjaz et al. 201 1 (Ic-BL), Levesque 
et al. 2010d (2009bb), and Levesque et al. 2010a (other GRB-SNe). Error bars 
illustrate measurement uncertainty, when published, plus a 0.07 dex systematic 
uncertainty. 

(A color version of this figure is available in the online journal.) 

The 39 LGRB host galaxies in the survey of Savaglio et al. 
(2009) are similarly low in mass and have high SFRs, with an 
average stellar mass of M* ~ 10 9 M 0 and SSFR ~ 3.5 Gyr -1 . 

The host galaxy of SN 2010ay falls below the L — Z relation 
for nearby star-forming galaxies, as illustrated by Figure 9. 
We have transformed the ^-corrected gri magnitudes from the 
MPA/JHU catalog to B band (Blanton & Roweis 2007). At 
the luminosity of the host galaxy of SN 2010ay, the median 
metallicity and standard deviation of the SDSS galaxies on the 
T04 scale are log(0/H) + 12 = 8.93 ± 0.17; the host galaxy of 
SN 2010ay falls in the 3rd percentile. In other words, the host 
galaxy of SN 2010ay is a 2 a outlier from the L — Z relation. 
Similarly, Levesque et al. (2010a) and Han et al. (2010) suggest 
that the host galaxies of LGRBs fall below the L — Z relation 
as defined by normal star-forming galaxies, BCGs, and the host 
galaxies of SNe Ic. 

Mannucci et al. (2011) have explained the offset of LGRB 
host galaxies from the M — Z relation as a preference for LGRB s 
to occur in host galaxies with high SFR, as characterized by 
the fundamental metallicity relation (FMR) of Mannucci et al. 
(2010; see also Lara-Lopez et al. 2010). Using the extended 
FMR for low-mass galaxies from Mannucci et al. (2011), 
the host galaxy of SN 2010ay should have a metallicity of 
log(0/H) +12 = 8.20 given its stellar mass and SFR. The 
FMR is calibrated to the Nagao et al. (2006) metallicity scale, 
which is similar to that of PP04 at this metallicity. Given the 
intrinsic scatter in the extended FMR on the order of ~0.05 dex, 
this value is consistent with the PP04 value we measure from 
our Gemini spectrum: log(0/H) +12 = 8.19. Kocevski & West 
(2011) similarly explain the offset of LGRB host galaxies from 
the M — Z relation as an SFR effect but suggest that the long 
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GRB host galaxies have even higher SFR than would be implied 
by the FMR. 

SN 2010ay is an example of an SN Ic-BL where the host 
galaxy is consistent with the M—Z relation for star-forming 
galaxies but deviates from the L — Z relation due to its low 
stellar-mass-to-light ratio (Section 6.1). Its 2 a discrepancy from 
the L — Z relation would be hard to explain as an SFR rate effect 
alone because among galaxies in the MPA/JHU catalog without 
AGNs (as defined above) and with Mb within 0.1 mag of the 
host galaxy of SN 2010ay, the host galaxy has an SFR in the 
26th percentile (<ler discrepancy). 

7. DISCUSSION 

SN 2010ay has all the hallmark features associated with 
GRB-SNe, and yet we find no evidence of a relativistic 
explosion to sensitive limits. We are able to place constraints 
on the energy, density, velocity, progenitor mass-loss rate, and 
gamma-ray flux of any GRB that may have been associated with 
it. In particular, we may rule out the association of a GRB that 
looks similar to any spectroscopically confirmed GRB-SN to 
date, except for XRF 060218. 

The low metallicity of the host environment of SN 2010ay 
may be suitable for GRB jet formation in the “collapsar” model, 
but our observations strongly constrain any relativistic outflow 
(Sections 4 and 5). In MacFadyen & Woosley (1999), a high 
rate of rotation in the core of the progenitor is required to power 
a relativistic jet. A low metallicity is prescribed to suppress 
the line-driven winds that would deprive the core of angular 
momentum. Apparently supporting this model, Stanek et al. 
(2006) found that the isotropic prompt energy release of the 
GRB-SNe decreases steeply with metallicity, and other surveys 
have found observational evidence for the preferential occur- 
rence of GRB-SNe in low-metallicity host galaxies (Fynbo 
et al. 2003; Prochaska et al. 2004; Sollerman et al. 2005; 
Modjaz et al. 2006; Wiersema et al. 2007; Christensen et al. 
2008; Modjaz et al. 2008; Levesque et al. 2010a; Chornock et al. 
2010; Starling et al. 2011). Challenging this view is the recent 
discovery of SN 2009bb, a broad-lined, engine-driven SN Ic 
found in a high-metallicity host environment (Soderberg et al. 
2010; Levesque et al. 2010d; Pignata et al. 2011). In SN 2010ay, 
we have found the opposite case — a broad-lined SN Ic found in 
a low-metallicity host environment, but without any indication 
(via either radio or gamma-ray emission) of a central engine. The 
existence of SNe 2009bb and 2010ay emphasizes that progeni- 
tor metallicity is not the key factor that distinguishes GRB-SNe 
from broad-lined SNe Ic without associated relativistic ejecta. 

We compare the absorption velocity of SNe Ic-BL and engine- 
driven SNe (GRB-SNe and SN 2009bb) to the metallicity 
of their host environments in Figure 10. This comparison 
emphasizes the diversity of explosion and host galaxy properties 
observed in both engine-driven SNe and SNe Ic-BL. The 
engine-driven SNe with the largest velocity gradients (2003dh, 
a = —0.9 ±0.1; 2009bb, a = —0.86 ± 0.08) occur at 
metallicities different by a factor of five (2003dh, log(0/H) + 
12 = 8.0; 2009bb, log(0/H) + 12 = 8.7 on the PP04 scale). 
Furthermore, the velocity gradient of engine-driven SNe seems 
to be uncorrelated with the velocity at late times; 2010bh 
( a = —0.22 ± 0.07) and 2006aj (a = —0.3 ± 0.2) have similar 
velocity gradients, but 2010bh had velocities ~ 1 0,000 km s' 1 
larger at 30 days after explosion. Among SNe Ic-BL from low- 
metallicity environments (log(0/H) + 12 < 8.5), there is a large 
range in both the characteristic velocity and velocity gradient 
04° = [10± 1, 9±2, 21 ±2] and a = [-0.5±0.2, -0.3±0.2, 
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Figure 10. SN absorption velocity, as traced by the Si n 16355 feature, vs. host 
galaxy oxygen abundance for SN 2010ay and other Ic-BL (red) and engine- 
driven explosions (blue) from the literature (as in Figure 5). The range of 
velocities hatched for each object comes from the velocity at 10 days and at 
30 days after explosion, according to the power-law fits presented in Figure 5. 
The oxygen abundance measurements using the PP04 03N2 diagnostic are from 
Levesque et al. (2010a; GRB-SNe), Sahu et al. (2009; SN2007ru), Young et al. 
(2010; SN 2007bg), and Modjaz et al. (2011; other SNe Ic-BL). The range 
of oxygen abundance hatched reflects the error bars quoted in the literature 
(when stated) plus the ~0.07 dex systematic error of the PP04 03N2 diagnostic 
(Kewley & Ellison 2008). 

(A color version of this figure is available in the online journal.) 

-0.4 ± 0.3] for SNe [2003jd, 2007bg, 2010ay]). For the three 
SNe Ic-BL at higher metallicities, the characteristic velocities 
tend to be lower and the velocity gradients tend to be stronger 
(v|f = [6±2,5±2, 10.3 ±0.7] anda = [-0.8 ±0.4, -1.5 ± 

0.7, -0.5 ± 0.1] for SNe [1997ef, 2002ap, 2007ru]). However, 
a larger sample is needed to exclude the possibility of SNe 
Ic-BL from super-solar metallicity environments that have high 
characteristic velocities or shallow velocity gradients. 

The fact that a GRB was not detected in association with 
SN 2010ay, despite its similarities to the GRB-SNe, could 
indicate that the relativistic jet produced by this explosion 
was “suffocated” before it emerged from the progenitor star 
(MacFadyen et al. 2001). In this scenario, the duration of the 
accretion event onto the newly formed central engine is short- 
lived and the jet post-breakout outflow is not ultrarelativistic. In 
the process of being suffocated, the jet transfers momentum 
to the ejecta such that the spectrum is broad lined and the 
ejecta velocity is very high, even at late times, as we observe 
(Section 3). However, the low host environment metallicity we 
measure for SN 201 Oay, which is similar to GRB-SNe, does not 
suggest high angular momentum loss that would help to weaken 
the jet. Another alternative is that binary interaction plays a 
key role in the commonality of relativistic outflows among 
SNe Ic-BL. 

Looking forward, additional SNe Ic-BL in sub-solar metal- 
licity host environments will be found through current and fu- 
ture generations of untargeted transient searches. Untargeted 
searches are required to find SNe from low-metallicity host en- 
vironments, because targeted searches only probe the highest 
metallicity galaxies due to the luminosity-metallicity relation- 
ship (Modjaz et al. 2011; Leloudas et al. 2011). 

8. CONCLUSIONS 

The optical photometric and spectroscopic, radio, and 
gamma-ray observations of SN 2010ay presented here pro- 


Table 5 

Comparison between SN 201 Oay and GRB 100316D/SN 2010bh 


Property 

SN 2010bh 

SN 201 Oay 

Host galaxy properties 

log(0/H)+12 a 

8.2 

8.19 

Redshift (z) 

0.059 

0.06717 

Mr 

-18.5 

-18.94 

Explosion properties 

SN type 

Ic-BL 

Ic-BL 

u 3 ° (10 3 km s -1 ) b 

24 ± 3 

21 ± 2 

Mr 

-18.60 ± 0.08 

-20.2 ± 0.2 

M Ni (M 0 ) 

0.10 ± 0.01 

0 9 +01 
u - y -o.i 

Me j (Mq) 

1.93-2.24 

>4.7 

E/c, 51 

12.0-13.9 

>10.8 

GRB energy (£ is0 . 

erg) >5.9 x 10 49c 

<6 x 10 48 


Notes. The observed properties of SN 2010bh and its host galaxy 
are given by Chomock et al. (2010), and light-curve modeling was 
performed by Cano et al. (201 la). The properties of SN 2010ay are 
derived in this paper. 

a The oxygen abundance derived from the PP04 03N2 metallicity 
diagnostic, as discussed in Section 6. 

b The absorption velocity at 30 days after explosion, as measured 
from the Si n 3.6355 feature in Section 3. 

c The lower limit of the total isotropic energy release estimated by 
Starling et al. (201 1). 


vide an example of a Type Ic-BL SN with explosion and host 
properties similar to the known GRB-SN 2010bh. This object 
demonstrates that SNe in low-metallicity environments with 
high-velocity ejecta are not necessarily accompanied by the tra- 
ditional signature of radio emission associated with long-lived 
relativistic jets. The existence of SN 201 Oay and SN 2009bb 
(a central-engine-driven event from a high-metallicity host en- 
vironment) indicates that progenitor metallicity may not be the 
key factor that distinguishes GRB-SNe from normal broad-lined 
SNe Ic. 

We conclude the following: 

1. Pre-discovery imaging of the SN 201 Oay from the Pan- 
STARRS1 3 Tt survey allows us to tightly constrain the 
early-time light curve of SN 201 Oay (see Figure 2) and 
explosion date (2010 February 21.3 ± 1.3), allowing us to 
search for gamma-ray emission that may have been asso- 
ciated with the explosion. By fitting the template SN Ib/c 
light curve of Drout et al. (2011), we derive an R-band 
peak absolute magnitude of —20.2 ± 0.2 mag — making 
SN 2010ay among the most luminous SNe Ib/c ever ob- 
served. This peak magnitude suggests that a large mass of 
nickel, /W N ; ~ 0.9* q 1 ; M 0 , has been synthesized. We esti- 
mate a ratio of My, to M c] that is ~2 times larger than in 
known GRB-SNe. 

2. Spectroscopy (see Figure 3) at the explosion site in the host 
galaxy of SN 201 Oay indicates that the host environment of 
the progenitor star had a significantly sub-solar metallicity 
(Z ~ 0.3 Z 0 ), similar to the host environments of known 
GRB-SN progenitors. 

3. The Type Ic-BL SN 201 Oay strongly resembles the 
GRB-SN 2010bh, particularly in light of its unusually high 
absorption velocities at late times (t>si ^ 19.2 x 10 3 km s' 1 
at 14 days after peak) and low-metallicity host environment. 
The comparison between these two SNe is summarized in 
Table 5. 
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4. Non-detections in late-time EVLA radio observations of 
the SN rule out the association of a GRB of the nature 
of the spectroscopically confirmed GRB-SNe, except for 
the radio afterglow associated with XRF 060218. Our 
radio observations imply limits on the velocity, energy, and 
density of any associated relativistic jet and the mass-loss 
rate of the progenitor (see Figures 4, 6-8). Additionally, no 
coincident gamma-ray emission was detected by satellites: 
the non-detection by the IPN indicates Ey< 6 x 10 48 erg, 
while the non-detection by the Swift-BAd indicates that the 
peak energy of the burst was <1 x 10 47 erg s -1 if the burst 
occurred during the ~20% of the explosion window when 
it was in the field of view of the instrument. This rules out 
associated prompt emission similar to that of GRBs 03 1203, 
030329, or 100316D, but not GRBs 980425 or 060218. 

The pre-discovery imaging of SN 2010ay demonstrates the 
capability of the untargeted PS1 survey for identifying and 
monitoring exotic transients, not only in its high-cadence 
medium-deep fields but also in the all-sky 37r survey. Additional 
detections and multi-wavelength follow-up observations of 
SNe Ic-BL will help to illuminate the role that ejecta velocity 
and progenitor metallicity play in the GRB-SN connection. 
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